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ABSTRACT 

We spectroscopically measure multiple hydrogen Balmer line profiles from laboratory plasmas to 
investigate the theoretical line profiles used in white dwarf atmosphere models. X-ray radiation 
produced at the Z Pulsed Power Facility at Sandia National Laboratories initiates plasma formation 
in a hydrogen-filled gas cell, replicating white dwarf photospheric conditions. Here we present time- 
resolved measurements of H/3 and fit this line using different theoretical line profiles to diagnose 
electron density, n e , and n = 2 level population, ri 2 - Aided by synthetic tests, we characterize the 
validity of our diagnostic method for this experimental platform. During a single experiment, we 
infer a continuous range of electron densities increasing from n e ~ 4 to ~ 30 x 10 16 cm” 3 throughout 
a 120-ns evolution of our plasma. Also, we observe ri 2 to be initially elevated with respect to local 
thermodynamic equilibrium (LTE); it then equilibrates within ~ 55 ns to become consistent with LTE. 

This supports our electron-temperature determination of T e ~ 1.3 eV (~ 15,000 K) after this time. 

At 7i e > 10 17 cm —3 , we find that computer-simulation-based line-profile calculations provide better 
fits (lower reduced X) than the line profiles currently used in the white dwarf astronomy community. 

The inferred conditions, however, are in good quantitative agreement. This work establishes an 
experimental foundation for the future investigation of relative shapes and strengths between different 
hydrogen Balmer lines. 

Keywords: line: profiles - plasmas - methods: laboratory: atomic - techniques: spectroscopic white 
dwarfs 


1. INTRODUCTION 

Recently, Tremblay & Bergeron (2009, hereafter TB) 
revamped white dwarf (WD) atmosphere models by 
modifying the constituent theoretical hydrogen line pro¬ 
files, which are calculated using the Unified Theory of 
Stark broadening (Smith et al. 1969; Vidal et al. 1970, 
1971, 1973, hereafter VCS). They followed the lead of 
Seaton (1990) to incorporate the occupation probability 
formalism of Hummer & Mihalas (1988) into the line- 
profile calculation to account for the loss of upper-energy- 
level transitions due to the high electric microfields com¬ 
mon in WD-atmosphere plasmas. This significantly im¬ 
proves the consistency of effective temperature, T e g, and 
surface gravity, log g, inferred from different, spectro¬ 
scopically observed hydrogen Balmer lines. It also results 
in a systematic increase in both T e g and log g. 

TB used their line profiles to spectroscopically deter¬ 
mine an increased mean mass of (M) = 0.649 Mq for 
DA WDs from the Palomar-Green Survey (Liebert et al. 
2005). Falcon et al. (2010a) used gravitational redshifts 
- a method mostly independent of theoretical line pro¬ 
files - to infer a DA mean mass of (M) = 0.647lg q) 3 Mq 
from the European Southern Observatory supernovae la 
progenitor survey (Napiwotzki et al. 2001). This is in 
excellent agreement with TB. Subsequent spectroscopic 
studies using TB line profiles, however, do not find mean 
masses from other WD samples that agree with that 
determined from Falcon et al. (2010a). These include 
Limoges & Bergeron (2010), Tremblay et al. (2011), and 
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Kleinman et al. (2013), who found ( M) = 0.606, 0.613, 
and 0.623 M 0 , respectively. 

The impact of theoretical line profiles on the interpre¬ 
tation of observed WD spectra as well as the discrepancy 
between the inferred mean masses from the spectroscopic 
and gravitational-redshift methods motivate the develop¬ 
ment of laboratory experiments capable of testing line- 
profile theories. TB profiles now replace the tabulated 
VCS profiles of Lernke (1997) as the WD-astronomy stan¬ 
dard. The difference between these calculations increases 
with increasing principal quantum number, n, and with 
electron density, n e . Therefore, measurements of the rel¬ 
ative line shapes of hydrogen Balmer lines formed at high 
electron densities can be used to discriminate between 
theoretical line-profile models. However, few laboratory 
experiments push into high n e (> 10 17 cm” 3 ) for the in¬ 
vestigation of H line shapes (e.g., Morris & Krey 1968; 
Parigger et al. 2003, 2008; Djurovic et al. 2005, 2009), 
and none of these measure multiple lines in the same 
plasma. Simultaneous measurements of multiple lines is 
highly desirable because it alleviates potential systematic 
uncertainties and promotes measurement accuracy. 

Our experiments create plasmas at WD photospheric 
conditions and reaches n e > 10 17 cm” 3 while simulta¬ 
neously observing multiple spectral lines (Falcon et al. 
2010b, 2013a, b; Montgomery et al. 2015). Hence, we 
provide the first such laboratory data at these plasma 
conditions; we do so using a fundamentally different ex¬ 
perimental approach than previous studies (e.g., Wiese 
et al. 1963; Hill & Gerardo 1967; Bengtson et al. 1969; 
Helbig & Nick 1981) by radiatively driving plasma forma¬ 
tion (e.g., Mancini et al. 2009; Liedahl 2011). Because we 
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use time-resolved spectroscopy, we also measure through¬ 
out a continuous range of n e in a single experiment and 
from the same plasma. 

Furthermore, we observe our plasma in absorption. 
This, along with simultaneously observing multiple lines, 
provides a unique opportunity to experimentally explore 
occupation probabilities (Hummer & Mihalas 1988) by 
measuring relative line strengths. 

In this paper we present spectroscopic measurements of 
the H/3 line. Historically, theoretical calculations, partic¬ 
ularly widths, of this line agree with measurements from 
benchmark laboratory experiments (e.g., Wiese et al. 
1972), validating its accuracy (for n e < 10 17 cm -3 ) to di¬ 
agnose plasma conditions (e.g., Kelleher et al. 1993). We 
thus fit H/3 to infer electron density, n e , and n = 2 level 
population, 712 , finding agreement between different the¬ 
oretical line-profile models (Section 3.3). These n e and 
7i2 measurements also allow us to infer electron temper¬ 
ature, T e (Section 3.4). We then scrutinize our sensitiv¬ 
ity to experimental uncertainties. Section 4.1 syntheti¬ 
cally investigates the dependence on T e , and Sections 4.2 
and 4.3 investigate inhomogeneities or gradients in the 
plasma conditions. Our scope is to secure our measured 
H/3 line as a diagnostic anchor for the investigation of 
relative line profiles (shapes and strengths) as part of the 
laboratory study of WD photospheres. 


source along a radial line of sight (LOS). X-rays irradiate 
the cell, propagate through a thin (1.44±0.02 /rm) Mylar 
window and through the H 2 gas, and are absorbed by a 
gold wall at the back end of the cell cavity. The gold 
heats to an electron temperature of T e ~ few eV and re¬ 
emits as a continuum, which photoionizes the hydrogen. 


CONTINUUM 



Figure 2. Cross-section false-color hardware drawing of the gas¬ 
cell central cavity. We illustrate the three observing lines of light 
(dashed, red lines). The Absorption LOS observes the same area on 
the back-lighting surface as does the Continuum LOS (vertical) but 
through the length of the plasma-filled cell. The Emission LOS also 
traverses this length but has no back-light. Thus, both horizontal 
lines of sight observe the same nominal plasma self-emission. 
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2. EXPERIMENTAL PLATFORM 

Ours is part of the Z Astrophysical Plasma Properties 
(ZAPP) Collaboration (Rochau et al. 2014), a group of 
experiments conducted simultaneously at the Z Pulsed 
Power Facility (e.g., McDaniel et al. 2002; Matzen et al. 
2005; Rose et al. 2010; Savage et al. 2011) at Sandia 
National Laboratories. ZAPP experiments take place in 
a large (> 60 m 3 ) vacuum chamber, each making use of 
the same 2 -pinch dynamic hohlraum (e.g., Spielman et al. 
1998; Nash et al. 1999; Stygar et al. 2001; Sanford et al. 
2002; Bailey et al. 2006; Rochau et al. 2008) x-ray source 
to initiate plasma formation. 



Figure 1. False-color hardware drawing of our “ACE”- 
configuration gas-cell assembly, capable of simultaneously observ¬ 
ing (1) the H plasma in Absorption, (2) the Continuum emis¬ 
sion from a gold, back-lighting surface, and (3) the H plasma in 
Emission. For one line of sight we show the hardware within the 
buffer (orange) and optics shield (gray). 

We place a gas-cell assembly (Figure 1), filled with hy¬ 
drogen (H 2 ), a distance away (324±2 mm) from the x-ray 


We observe the H plasma along lines of sight parallel 
to the gold wall and perpendicular to the photoioniz- 
ing radiation (Figure 2). For each LOS, a 50-mm-long 
buffer cavity separates the optics from the plasma formed 
within the central cavity. Optical fibers deliver the light 
to time-resolved spectrometer systems. Sections 4.2 and 
4.3 elaborate on our observations (Figures 12 and 14 il¬ 
lustrate). See Falcon et al. (2013b) for a description of 
our experimental platform and plasma formation, and 
see Falcon et al. (2015) for details of our data processing 
and calibration. 

3. TIME-RESOLVED SPECTROSCOPIC 
MEASUREMENTS 

Inside the gas-cell central cavity a polyhedral stainless 
steel block, coated with 5 /im of gold, rests on one end. 
Its surface is tilted with respect to the plane normal to 
the 2 -pinch x-rays ( y-z ) and pitched with respect to the 
horizontal {x-y). This allows the x-rays, the optics in 
one horizontal LOS (Absorption), and the optics in the 
vertical LOS (Continuum) to each have a direct view. 
The x-rays heat this surface like they do the gold back 
wall, allowing it to serve as a back-light for absorption 
measurements. 

3.1. Extracting Transmission 

Absorption spectra of plasmas contain three compo¬ 
nents: 

jmeas jback—light ^plasma ^plasma /-j\ 

/back-light is the spectral radiance of the back-lighting 
continuum. The transmission, j , P lasma j anc j the self¬ 
emission, jP lasma ^ both describe the plasma but in differ¬ 
ent ways. We focus on yP lasma f or our analysis because 
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of its preferable signal-to-noise (S/N) ratio. It also pro¬ 
vides a constraint on relative line strengths because all 
the transitions of the absorption lines we observe share 
the same initial (lower) state (n = 2); this is a subject of 
future work. 

The components in Equation 1 are separately mea¬ 
surable using the “ACE” gas cell. Because our absorp¬ 
tion LOS ( I x bs ) observes a back-lighting surface through 
a length of plasma, I x hs equals /™ eas . Our contin¬ 
uum LOS, which observes this same back-lighting sur¬ 
face (Figure 2), also measures absorption spectra. How¬ 
ever, since the length of intervening plasma is minimal 
(L ~ 7 mm), the transmission is close to unity (absorp¬ 
tion is small) and we assume negligible self-emission, so 
that I c x ° nt »i£ ack - light . 

With our emission LOS we approximate I x m to equal 
^plasma i en gth 0 f plasma this LOS observes (L ~ 

120 mm) nearly matches that of the absorption LOS {L ~ 
114 mm), which is ~ 5 % shorter. We also approximate 
the nature of the plasma emission from each LOS to be 
the same, though the absorption LOS includes a ~ few- 
mm region adjacent to a gold, back-lighting surface that 
the emission LOS does not (Figure 2). 

Upon placing our measured emission and absorption 
data onto the same absolute scale (Falcon et al. 2015), 
we correct the latter for plasma self-emission. We can 
then substitute our measurements into Equation 1 and 
invert to extract transmission: 

Tabs,cor Tabs _ rem 

^plasma ^ jA _ _A_/^j 

± A r '" / rcont ~ rcont ’ v ' 

1 \ 1 \ 

Figure 3 shows an example H/3 corrected-absorption 
spectrum (solid, red; j abs ’ cor ) versus its uncorrected spec¬ 
trum (dotted, red; I x hs ) for an integration of the time- 
resolved data over a 10-ns interval (line-out). 


z2553, 55 ns after x-ray onset 



Wavelength (A) 

Figure 3. Example H j3 corrected-absorption spectrum (solid, red) 
plotted with its uncorrected-absorption (dashed, red) and with the 
subtracted-emission (dashed, blue). These 10-ns line-outs center 
at 55 ns after the onset of x-rays from experiment z2553. We 
also show the back-lighting continuum shape measured from z2554 
and normalized to the z2553 corrected-absorption. This contin¬ 
uum shape (dotted, green) is similar to the straight-line continuum 
(solid, green) we use to convert these data to transmission. 

A desirable strategy to address plasma self-emission 
for absorption data measured for opacity studies is to 


require a back-lighter that is sufficiently bright to over¬ 
whelm self-emission (e.g., Davidson et al. 1988; Perry 
et al. 1996; Bailey et al. 2007, 2009, 2015). Our simulta¬ 
neous emission and absorption measurements shed this 
requirement by allowing us to remove the former from the 
latter, and since we are the first to study H line profiles in 
absorption, we are also the first to utilize this technique 
for H. Consequently, this strategy, and in particular the 
approximations we list, deserve further scrutiny. 

At the electron densities we encounter, overlap of the 
H/3 wings with adjacent lines (i.e., Ha and Hy) is mini¬ 
mal. Also, we do not observe the bound-free continuum 
(e.g., Dappen et al. 1987) to reach far enough redward 
to affect H/3. To focus on the H/3 bound-bound transi¬ 
tion we simplify our transmission extraction by redefin¬ 
ing the continuum baseline (/£ ont ). Since our measured 
shapes of the back-lighting continua (using the contin¬ 
uum LOS 3 ) are quite linear in the H/3 spectral vicin¬ 
ity (dotted, green in Figure 3), we use straight-line fits 
(solid, green); /™ nt ss 7“ ont . This circumvents potential 
contamination due to any additional continuous (weakly 
dependent on wavelength) opacity sources, such as that 
due to the H _ ion (e.g., Wishart 1979; Griem 1997). 

3.2. Spectral-Line Fitting Model 

Our fitting model is the transmission due to bound- 
bound absorption of a hydrogen Balmer line for a homo¬ 
geneous (single T e , n e , and n\) plasma: 

T x~ h = ,t P*exp{-K h x ~ h L} , (3) 

where if) is the measured instrumental broadening, * is 
the convolution operator, K b_b is the opacity due to pho¬ 
toexcitation for a bound-bound transition while neglect¬ 
ing stimulated emission (e.g., Mihalas 1978), and L is the 
length of the observed plasma. 

We express the opacity as 

7re 2 

K b ~ b =- f\^. u niw u (n e )<t)\(T e ,n e ), (4) 

m e c 

where e is the electron charge, m e is the electron mass, 
c is the speed of light, fi-> u is the oscillator strength of 
the transition from the lower to upper energy level, and 
n\ is the lower-level population. The upper-level occu¬ 
pation probability, w u (n e ), depends on electron density 
(Hummer & Mihalas 1988). The line profile, <j)\(T e ,n e ), 
is normalized (f 0\dA = 1) and a function of both elec¬ 
tron temperature, T e , and n e . Because the dependence 
on T e is relatively weak (Section 4.1) and because we 
have empirical constraints (Section 3.4), our fits assume 
T e = 1.0 eV. 

Combining the constants into a single factor, C K , and 
specifying the H/3 transition, Equation 3 becomes 

T x ^ = ip* exp ^-C K f 2 ^. 4 n 2 w 4 (n e )(/) x /3 (n e )Lj . (5) 

3 Because only two spectrometer systems are sufficiently cali¬ 
brated (Falcon et al. 2015), we do not use the continuum LOS 
to measure /J ont ^/back-light^ Qn ex p er i men t z 2553. This experi¬ 
ment uses one spectrometer system for the absorption LOS and the 
second for the emission LOS, enabling the self-emission correction. 
On the subsequent experiment, z2554, the second spectrometer sys¬ 
tem uses the continuum LOS to measure the back-lighting spectral 
shape needed to infer transmission. 
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The principal quantum numbers of the initial and final 
states are n = 2 and n = 4, respectively. 

We use oscillator strengths from Baker (2008), who 
reports to high precision (< 0.3% uncertainty). The oc¬ 
cupation probability comes from Seaton (1990); for H/3, 
w u (n e ) is quite nearly unity. We fix L according to the 
dimensions of our gas-cell hardware design. This leaves 
ri 2 and cj^f (n e ) to be determined from experiment. 

3.3. Fits to Measured H/3 Transmission 

Figures 4 and 5 plot twelve consecutive 10-ns line-outs 
of H/3 transmission (black diamonds) we measure from 
experiment z2553; the first line-out begins at the onset 
of x-rays (centered at 5 ns after onset). Uncertainties in 
spectral points (black, vertical lines) reflect the S/N. 

We use a Levenberg-Marquardt minimization (Leven- 
berg 1944; Marquardt 1963) to fit our measured trans¬ 
mission spectra with Equation 5. Red, blue, pink, and 
green curves correspond to fits using theoretical line pro¬ 
files, , that follow Vidal et al. (1973, VCS), Tremblay 
& Bergeron (2009, TB), Gigosos et al. (2003, GGC), and 
Gomez et al. (2015, Xenomorph or XENO), respectively. 
The VCS profiles come from Lemke (1997). 

These theoretical line profiles are quite similar, differ¬ 
ing noticeably only in the line center (enlarged window 
in lower left-hand side of each panel). Fits using VCS 
and TB both based on the Unified Theory - are nearly 
identical. Fits using GGC and Xenomorph - both incor¬ 
porating computer simulations (e.g., Stamm et al. 1984; 
Stambulchik & Maron 2010) - exhibit less-pronounced 
central humps 4 , and this latter theory is the only one of 
the four that calculates asymmetry. 

We indeed measure structure in the line center, includ¬ 
ing asymmetry (e.g., Kudrin & Sholin 1963; Wiese et al. 
1975; Halenka 1988; Djurovic et al. 2005). This struc¬ 
ture remains absent at the lowest n e and increases with 
this parameter. Early and late in the experiment the fits 
reproduce our measurements near the line center quite 
well, but from 35-95 ns there are statistically significant 
deviations. The fits predict central structure that we do 
not observe over this duration of the experiment, which 
spans the range of 1 x 10 17 < n e < 3 x 10 17 cm -3 . Fur¬ 
ther work is required to ascertain whether this is due 
to deficiencies in the theoretical line profiles or in the 
experiment. 

GGC and Xenomorph profiles provide systematically 
better fits, as expressed by reduced Xred (Figure 6). 
Though the VCS and TB profiles yield poorer fits, all 
theories infer values for electron density, n e , and lower 
(n = 2) level population, 712 , that agree (Figures 7 and 
8, respectively). The standard deviation of n e between 
theories increases from < 2 to ~ 6 % of (n e ) as time in¬ 
creases (as n e increases). For n%, the trend is the same, 
increasing from < 2 to ~ 4 %. 

We monitor n e increase throughout the first 95 ns of 
experiment z2553 before reaching a plateau of n e ~ 
3 x 10 17 cm~ 3 . Beyond 55ns we do not plot n e inferred 
using GGC profiles, because these values extend beyond 
the n e range covered by Gigosos et al. (2003). In the 
first 25 ns we overlap (4 x 10 16 < n e < 10 17 cm -3 ) the 

4 The central feature of a line profile is usually described as the 
central dip because it is usually observed in emission. In our novel 
approach, however, we observe the profile in absorption. 


electron-density range observed by Wiese et al. (1972) - 
the only other study to measure multiple H Balmer lines 
near these conditions. We then exceed it by approxi¬ 
mately a factor of three. 

The right-hand-side, vertical axis of Figure 7 shows 
the ionization, which we define as —?¥—,where n° is 
the neutral H atom density. To determine the total par¬ 
ticle density (ntot = n e + n°), we measure the initial 
H 2 -gas fill pressure inside our cell in situ using piezore- 
sistive pressure sensors (see Appendix A in Falcon 2014). 
Assuming that the photoionizing radiation forming our 
plasma completely dissociates the H 2 , the Ideal Gas Law 
translates our measured pressure (P = 10.63 ± 0.05 Torr) 
to ntot = (6.92 ± 0.08) x 10 17 cm -3 . Uncertainties for P 
and ntot reflect the precision of the pressure-sensor volt¬ 
age measurement. They do not include all systematic 
uncertainties, such as that due to the accuracy of the 
correction for zero-pressure output voltage. 


3.4. Electron-Temperature Determination 

The lower-level population ( 712 ) also increases with 
time after exhibiting a moderate spike soon (~ 20 ns) af¬ 
ter the onset of x-rays (Figure 8). To give our measured 
712 context, we compare it to that of a plasma in local 
thermodynamic equilibrium (LTE) by plotting the ra¬ 
tio (open diamonds connected by dotted lines and corre¬ 
sponding to the purple, right-hand axis). We determine 
the LTE n = 2 population by first inserting our inferred 
electron density and measured total particle density into 
the Saha equation and solving for electron temperature. 
We neglect ionization potential depression (e.g., Crowley 
2014) and approximate all neutral atoms to be in the 
ground state (n° s=s rii): 



T^tot ^1 


/ 2nm e kT e \ 3 ^~ 

) 6XP 



( 6 ) 


where k is the Boltzmann constant, h is the Planck 
constant, and the neutral-hydrogen ionization energy is 
e° = 13.6 eV. Then, using this T e with the Boltzmann 
relation yields ti 2 (LTE): 


7 i 2 (LTE) f -(e 2 - ei) 1 

-= —exp -—-> , 

m g 1 ( kT e J 


where the statistical weight, g n , equals 2n 2 , and the en¬ 
ergy difference between the n = 2 and 71 = 1 states is 
e 2 — ei = 10.2 eV. Early in time, our measured 712 is ~ 5 
times that at LTE (dashed, purple, horizontal line). This 

ratio ( ra2 ^™(LTE) ed ^ ) monotonically approaches unity and 
reaches it ~ 55 ns after the onset of x-rays. This implies 
that the ionization fraction, n = 2 population, and n = 1 
population are consistent with their LTE values. 

Once in LTE, T e ~ 1.25 eV (~ 14,500 K) and rises 
to T e ~ 1.35 eV (~ 15, 700 K) by ~ 115 ns (Figure 9). 
Before ~ 55 ns, our LTE T e estimate provides an upper 
limit (open diamonds, dotted lines) to the true T e , as¬ 
suming our plasma is “overionized” (e.g., Kawasaki et al. 
2002). Our measurement of “overexcitation” (i.e., ele¬ 
vated n = 2 population) supports this assumption, as 
does the premise that our plasma is initially photoion- 
ized (Falcon et al. 2013b). 
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5 ns z2553 15 ns 



45 ns 55 ns 




Wavelength (A) 

Figure 4. Time sequence of H/3 transmission spectra (black diamonds) measured from our laboratory plasma during experiment z2553. 
Uncertainties (black, vertical lines) reflect the signal-to-noise ratio. The first of these 10-ns line-outs begins at the onset of x-rays to the gas 
cell. Fits (solid curves) use the theoretical line profiles of Vidal et al. (1973, VCS; red), Tremblay Sz Bergeron (2009, TB; blue), Gigosos 
et al. (2003, GGC; pink), and Gomez et al. (2015, XENO; green). We look closely at the line center where the fits differ most. 
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65 ns Z2553 75 ns 



105 ns 115 ns 




Wavelength (A) 


Figure 5. Same as Figure 4 but plotting the subsequent six 10-ns line-outs. Here we do not include GGC fits since the inferred electron 
densities exceed the range spanned by the GGC line-profile grid. 


sources, which we add in quadrature. cr fit is from the 
fit, which is a random uncertainty due to noise in the 
spectra. It is typically ~ 2 % and ~ 1 % for n e and 
ri 2 , respectively. cr cal is due to the transmission extrac- 


3.5. Uncertainties in Inferred Plasma Conditions 

The uncertainties for our inferred n e (Figure 7), 712 
(Figure 8), and, subsequently, T e (Figure 9) have two 
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z2553 



Figure 6. Goodness-of-fit (reduced chi squared, X^ed) versus time 
for fits to H/3 transmission measured from experiment z2553 and 
using VCS (red), TB (blue), GGC (pink), and Xenomorph (green) 
line profiles. Early in time when n e is low, all theories provide 
equally good fits. Later in time when n e > 10 1 ' cm -3 , GGC and 
Xenomorph profiles provide better fits than VCS and TB. 


z2553 



Time (ns) 

Figure 7. Electron density (left-hand-side, vertical axis) and ion¬ 
ization (right-hand-side, vertical axis) versus time determined from 
fits to measured H/3 transmission. Section 3.5 describes uncertain¬ 
ties (vertical lines). We time-resolve a smooth increase in electron 
density (n e ) which allows us to measure line profiles throughout a 
range of n e in a single experiment and from the same plasma. All 
theories infer n e in agreement. 

tion (Section 3.1), which depends upon the accuracy of 
our calibration (i.e., the relative measurement of plasma 
emission and absorption from different spectrometer sys¬ 
tems; Falcon et al. 2015). 

Early in time, the magnitude of <r cal is similar to that of 
a Rt . Throughout the experiment, however, a cal increases 
to ~ several %. This is a result of the evolution of the 
back-lighting continuum emission (J™ nt ) which decreases 
as the gold surface cools, thus increasing the significance 
of the self-emission component (J| m ) to the measured 
absorption (/ abs ). 

3.6. Plasma Reproducibility 

We observe the same qualitative trends in goodness- 
of-fit, electron density (Figure 10), and lower-level pop¬ 
ulation versus time across multiple experiments, giving 
credence to plasma reproducibility using our experimen¬ 


z2553 



Figure 8. Left-hand-side axis: similar to Figure 7 but for lower 
(n = 2) level population versus time. Right-hand-side axis: ratio 
of measured versus LTE n = 2 population (open diamonds, dotted 
lines). By ~ 55 ns after the onset of x-rays the n = 2 population is 
at its LTE value (dashed, purple, horizontal line). 

z2553 
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Figure 9. Electron temperature ( T e ) versus time inferred from the 
Saha equation assuming local thermodynamic equilibrium (LTE) 
and using our measured ionization fraction (diamonds). Before the 
plasma reaches LTE at ~ 55 ns (Figure 8), these inferred T e values 
are upper limits assuming the plasma is “overionized” (see text). 

tal platform. This includes the ri 2 spike at ~ 20 ns. 

Experiments z2554 (blue) and z2588 (green) have no 
/| m measurement, so we infer n e from fits to 7j( bs (i.e., 
transmission uncorrected for plasma self-emission; open 
diamonds, dashed lines). For z2553 (red), foregoing the 
emission correction overestimates n e by ~ few % early 
in time, increasing to ~ few tens of % late in time. It 
also underestimates ri 2 by percentages similar to the n e 
overestimations. 

4. SYNTHETIC INVESTIGATIONS INTO 
SYSTEMATIC UNCERTAINTIES 

Our fitting procedure (Section 3.2) approximates the 

line profile, , to have a negligible dependence on elec¬ 
tron temperature, T e , and therefore depend solely on 
electron density, n e . It also assumes that our experi¬ 
mental plasma either (1) is homogeneous and exists at 
a single plasma condition or (2) can be sufficiently and 
uniquely described as a homogeneous plasma. 
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Figure 10. Similar to Figure 7. The inferred electron density ver¬ 
sus time from multiple experiments is qualitatively reproducible. 
Also, fits to H/3 transmission uncorrected (open diamonds, dashed 
lines) for plasma self-emission overestimate the electron density 
compared to fits with the correction (filled diamonds, solid lines). 
We show fits using VCS theoretical line profiles. 

We implement synthetic data analysis (e.g., Nagayama 
et al. 2012a, b) to investigate the sensitivity of our fits to 
the aforementioned factors. We simulate H/3 transmis¬ 
sion data at different T e (Section 4.1) as well as those 
resulting from inhomogeneous conditions (Sections 4.2 
and 4.3). These data use VCS theoretical line profiles, 
and include random noise (S/N = 33) and instrumen¬ 
tal broadening (Voigt profile with ~ 10-A full-width-at- 
half-maximum) typical for our experiments. They also 
span a range of plasma conditions (i.e., 10 16 < n e < 
4 x 10 17 cm -3 in steps of A n e = 3 x 10 16 cm -3 and 
3 x 10 14 < n? < 1.2 x 10 15 cm -3 in steps of Ari 2 = 10 14 
or 3 x 10 14 cm -3 ). We independently adjust T e , n e , and 
7i2 without assuming LTE. 

4.1. Sensitivity to Electron Temperature 

Our ionization fraction (-^) and n -2 measurements 

indicate that our plasma reaches LTE by ~ 55 ns (Fig¬ 
ure 8). This corresponds to a thermalized temperature 
of T e ~ 1.3 eV. Before this time and while our plasma is 
not in LTE, we assume that our T e estimates are upper 
limits (Figure 9); a precise T e determination is nontrivial, 
likely requires collisional-radiative hydrodynamic model¬ 
ing (e.g., Hansen et al. 2007), and is beyond the scope of 
this paper. Because of this and because of the relative in¬ 
sensitivity we show in Figure 11, we adopt a compromise 
(as a function of time and hence T e ) by fixing T e = 1.0 eV 
in our spectral-line fitting model (Section 3.2). 

By assuming T e = 1.0 eV, our fits underestimate n e 
(top panel) and ri 2 (bottom panel) for cooler ( T e < 
1.0 eV) plasmas and overestimate for warmer (T e > 
1.0 eV) plasmas. This effect is more significant for n e 
than for 712 , for which it is quite small (< 1 % for most 
cases). Uncertainties (vertical lines) only simulate er fit 
and do not include er cal (Section 3.5). 

4.2. Gradient in Electron Density Transverse to our 
Line of Sight 

As a consequence of the plasma formation inside our 
gas cell, ionization decreases with increasing distance 




Figure 11. Ratio (vertical axis) of inferred parameter (fit) and the 
value of that parameter for synthetic H/3 transmission (synth) cre¬ 
ated at different electron temperatures (horizontal axis). Our fits 
assume T e = 1.0 eV. Top and bottom panels show inferred electron 
density, n e , and lower (n = 2) level population, n 2 , respectively. 


from the gold wall (Falcon et al. 2013b). Since we ob¬ 
serve our plasma parallel to the gold wall (perpendicular 
to the photoionizing radiation) along lines of sight with 
finite diameter (~ 3 mm; Falcon et al. 2015), our mea¬ 
sured spectra sample some range of plasma conditions. 


Gold Wall 


Z 

♦ 
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Side View 



Figure 12. Cross-section schematic of the gas-cell central cavity 
from the point of view of a horizontal line of sight. The collection 
beam (red) traverses the hydrogen plasma parallel to the gold wall 
and perpendicular to the photoionizing radiation. We construct 
synthetic data simulating a transverse gradient in electron density, 
n e , by assembling the collection beam out of zones with different n e 
linearly decreasing with distance from the gold wall (x direction). 


We simulate a linear gradient in n e transverse (x di¬ 
rection) to our LOS using the mean of synthetic homo¬ 
geneous transmission spectra of different n e weighted ac¬ 
cording to fractional areas of a circle. This corresponds 
to our cylindrical collection beam. We describe the gra¬ 
dient as the difference of n e at the boundaries of the LOS 
collection beam divided by the mean n e : 

Gradient s K(i = "0 - -.(> = 1)1 x 100 %> (8) 

if E »,(•) 

i= 1 
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where zone i is one of m zones (Figure 12). 




Figure 13. Vertical axis: ratio of inferred parameter (fit) and 
the value of that parameter used to create synthetic H/3 transmis¬ 
sion (synth) that simulates a linear gradient in electron density, 
n e , transverse (x direction) to the observed line of sight (LOS). 
Horizontal axis: n e gradient transverse to synthetic LOS defined 
in Equation 8. Top and bottom panels show inferred n e and lower 
(n = 2) level population, 712 , respectively. The steeper the gradi¬ 
ent, the more the fit underestimates each parameter. 


The top and bottom panels of Figure 13 plot our results 
for n e and 712 , respectively. We show the ratio (colored 
diamonds) of the parameter inferred from fitting and the 

synthetic value the fit sought to recover (i.e., n "synth) 

and n ”( 2 S ynth) ) • For n e , this synthetic value is the mean 
n e across zones, and, since the gradient is linear, it equals 
that of the central zone. 

As the gradient increases, fits underestimate both n e 
and 712- Empirically, we find that neither nor 

^"(syntii) depend on n e , so we average over this param¬ 
eter. They do depend on 712 , however. As 712 (line 
strength) increases, fits underestimate n e less and 712 
more. 

Though an exponential function better describes the 
decrease in measured n e with increasing distance from 
the gold wall (Falcon et al. 2013b), we approximate this 
gradient to be linear across our 3-mm-diameter collection 
beam. The measured gradient from experiments z2300 
and z2302 is then ~ 100 % for a LOS centered at 10-mm 
away from the gold wall and using a ~ 30-Torr H 2 gas fill 
(Falcon et al. 2013b). Our synthetic test suggests that 
this could lead to an underestimate of both n e and 712 of 
up to ~ 3 %. 

Experiment z2553 uses a ~ 10-Torr gas fill. Lower 
pressure results in decreased attenuation of the photoion- 
izing radiation with increasing distance from the gold 
wall. This means we encounter higher n e for a fixed LOS 
distance. It also means the transverse gradient is not as 
steep. Therefore, the measured ~ 100-% gradient should 
be an overestimate for the z2553 observing LOS (dotted, 
purple, vertical line). The underestimation of n e and 712 
is then < 3 %. 


An inevitable attribute of laboratory plasmas is their 
finite size. Boundary layers exist where conditions tran¬ 
sition away from that of the bulk plasma, contribut¬ 
ing to inhomogeneity along the LOS. Observing long 
plasmas minimizes this effect (e.g., Bengtson & Chester 
1976), and though we create the longest photoionized 
(radiation-driven) laboratory plasmas to date (Falcon 
et al. 2013b), their lengths remain finite. 


Top View 


x 
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Figure 14. Top-view cross-section schematic of the gas-cell cen¬ 
tral cavity. We construct synthetic data simulating a line-of-sight 
gradient in electron density by appending a foreground boundary 
layer (diagonally striped, red) to the homogeneous-plasma compo¬ 
nent (solid, red) within a collection-beam volume. 


We assume that the boundary-layer plasma extending 
into the buffer cavity (Figure 1) is the dominant source 
of inhomogeneity along our observing LOS (y direction). 
Thus, we create synthetic data consisting of two compo¬ 
nents: a homogeneous plasma (solid, red in Figure 14) 
the length (L ~ 114 mm) of the absorption LOS inside 
our gas-cell central cavity, and a boundary layer (diag¬ 
onally striped, red) with variable length whose n e lin¬ 
early decreases along that length toward the observer. 
This homogeneous component neglects the region adja¬ 
cent to the gold, back-lighting surface (Section 3.1), and 
the boundary-layer n e does not decrease to zero from the 
electron density, n e (synth), of the homogeneous compo¬ 
nent. It decreases to the lowest value of our theoretical 
line-profile grid, n e = 10 15 cm -3 . 




Boundary Layer (mm) 

Figure 15. Vertical axis: ratio of inferred parameter (fit) and the 
value of that parameter used to create synthetic H (3 transmission 
(synth) that simulates a boundary-layer plasma in the foreground 
of a homogeneous plasma. Horizontal axis: length of synthetic 
boundary layer. The longer this layer, the more each fit underesti¬ 
mates n e (top panel) and overestimates 712 (bottom panel). 


4.3. Gradient in Electron Density Along our Line of . , , , , . r , r . 

Sight or a Boundary-Layer Plasma A s y nthetl <; boundary layer m the foreground of a ho¬ 

mogeneous plasma causes fits to underestimate n e and 
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overestimate 712 (Figure 15). Longer boundary layers ex¬ 
acerbate these effects. Empirically, they do not depend 
on ri 2 (similarly to the lack of n e -dependence of n ”sy„*h) 

and n^( 2 S ynth) f° r the transverse gradient in Section 4.2), 
so we average over this parameter. As n e (synth) in¬ 
creases, fits underestimate n e more and overestimate ri 2 
less. 

With current measurements, the true conditions in the 
buffer cavity are difficult to ascertain. As a conserva¬ 
tive estimate, though, we do not expect a boundary-layer 
plasma to extend beyond ~ 15 mm (dotted, purple, ver¬ 
tical line) into the 50-mm-long buffer. In similar exper¬ 
iments we measure n e to significantly fall off at 15 mm 
from the gold wall (Falcon et al. 2013b). This is in the 
x direction, perpendicular to the gold wall, where the 
LOS to the photoionizing radiation is direct, thus maxi¬ 
mizing the irradiance. Not only is the LOS to a bound¬ 
ary layer peripheral, a 5-mm-diameter aperture parti¬ 
tions the central cavity from the buffer cavity (Figures 2 
and 14), choking the radiation that can penetrate. At 
15 mm our synthetic test suggests that we underestimate 
n e and overestimate 712 by < 14%. 

5. CONCLUSIONS 

We examine hydrogen Balmer-/3 spectral line profiles 
measured from laboratory plasmas at white dwarf photo- 
spheric conditions; we establish this line as a diagnostic 
anchor for our experiments. This sets the stage for the 
analysis of the relative line shapes and strengths of multi¬ 
ple H Balmer lines we measure simultaneously (H/3, Hy, 
and H(5). 

The conditions we encounter during a single experi¬ 
ment smoothly increase throughout a 120-ns plasma evo¬ 
lution and span from n e ~ 4 to ~ 30 x 10 16 cm -3 , ex¬ 
ceeding, by a factor of three, the highest n e achieved 
in the benchmark work of Wiese et al. (1972). At 
n e > 10 17 cm -3 , we find that the theoretical line pro¬ 
files adopted throughout the WD astronomy community 
(i.e., Lemke 1997; Tremblay & Bergeron 2009) do not fit 
our measured H/3 profiles as well as computer-simulation- 
based calculations (i.e., Gigosos et al. 2003; Gomez et al. 
2015). Despite this, all these profiles infer similar n e and 
n = 2 level populations, 712 ; the standard deviations of 
inferred parameters between theories remain modest, in¬ 
creasing, as n e increases, from < 2 to ~ 6 % of (n e ) and 
from 2 to ~ 4 % of ( 712 ). 

We investigate the sensitivity of our spectral-line fit¬ 
ting model to the following factors: (1) the approxima¬ 
tion that the H/3 line shape has a negligible dependence 
on electron temperature, (2) a gradient in n e transverse 
to our observing line of sight (LOS), and (3) an interlop¬ 
ing boundary-layer plasma along our LOS. 

For each of these we have experimental guidance. (1): 
By monitoring 712 we witness our photoionized plasma 
smoothly relax into local thermodynamic equilibrium 
(LTE) in ~ 55 ns, which corresponds to T e ~ 1.3 eV 
(~ 15,000K). Since our fits assume T e = 1.0eV, we may 
overestimate n e by < 10 % (at our highest n e ) and 712 by 
<3%. (2) and (3): Using this experimental platform, 
we measure n e of a similar H plasma along lines of sight 
at different distances from the radiating gold wall (Fal¬ 
con et al. 2013b). This constrains the transverse gradient 


across our LOS for this experiment to be < 100 %. It also 
suggests that a boundary-layer plasma along our LOS 
is, conservatively, < 10 % of our homogeneous-plasma 
length. Due to (2), we should underestimate both n e 
and ri 2 by < 3%. Due to (3) we should underestimate 
n e and overestimate 712 by < 14%. 

Our ability to measure multiple H lines simultaneously 
at high electron densities (n e > 10 17 cm -3 ) provides us 
with a sensitive and unprecedented test of theoretical 
line profiles. Relative H Balmer line shapes can discrim¬ 
inate between theories (e.g., Vidal et al. 1973; Tremblay 
& Bergeron 2009); because we measure in absorption , 
relative line strengths will offer a new perspective on oc¬ 
cupation probabilities (Hummer & Mihalas 1988). 
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